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SUMMARY 

The differential equation governing the deflection of a disk of var- 
iable thickness subjected to gyroscopic loading is derived. For the case 
of a disk of constant thickness, solutions are obtained by a finite- 
difference method for a range of centrifugal loading parameter M from. 

0 to 50, and the ratio of shaft to disk radii varying from 0.1 to 0.3. 
Results are presented in dime nsionless form suitable for design purposes. 
The method for solving the problem of a disk of variable thickness with 
a temperature gradient is also presented. 


INTRODUCTION 

When the axis of rotation of a disk is itself rotated, forces are 
set up normal to the disk. These forces, or gyroscopic loads, occur on 
jet-engine compressor and turbine disks whenever the airplane changes 
direction either in the air or on the ground. These gyroscopic forces 
will deflect a disk out of its plane of rotation, induce vibratory bend- 
ing stresses, and produce a bending moment that will increase the shaft 
bearing loads . This problem may be even more serious for the high-speed 
disks used in missiles undergoing high accelerations. Although the 
stresses thus produced may not themselves be very large, the combined 
effect of these stresses with the already existing stress distribution 
may be sufficient in scsne cases to cause a failure . It is therefore nec- 
essary to investigate both the stresses and the deflections of a disk 
under gyroscopic loading. 

By considering the gyroscopic forces and centrifugal forces acting 
on an element of a disk of variable thickness, the general differential 
equation describing the deflection is derived. A finite-difference solu- 
tion of this differential equation is obtained for a constant -thickness 
disk rotating about its center at a constant angular velocity. The com- 
puted deflections and the corresponding stresses are plotted in dimen- 
sionless form for a wide range of a centrifugal- loading parameter, and 
for a number of different ratios of shaft to disk radii. 
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PROCEDURE 

The general differential equation governing the deflection of a 
disk of variable thickness subjected to gyroscopic loading is derived in 
appendix B and the solution of this general equation is discussed in 
RESULTS AND DISCUSSION . In this section the procedure used for obtain- 
ing numerical solutions of this equation is presented for the speci al 
case of a constant-thickness disk rotating about its center at a constant 
angular velocity cd and at the same time about a diameter fixed in space 
at an angular velocity ft and an angular acceleration Cl. (The symbols 
used in this report are defined in appendix A. ) Figure 1 shows such a 
disk along with the type of loading that results from such rotations . 

If the deflection at any point such as P is given by w for this 
disk, and if the dimensionless deflection for a given angle 9 is de- 
fined by 


Y 


w 

C 


where 

C = 2mn a 5 sin 0 + cos e) (l) 

and m, a, h, and D are the mass density, rim radius, disk thickness, 
and flexural rigidity, respectively, and 9 is the angular coordinate 
on the face of the disk measured from the diameter about which the disk 
is rotating, the general differential equation (B19) describing the 
dimensionless deflection Y becomes 



- ~ +M(l - p^ 

Lp> 



p 


( 2 ) 


where p is the dimensionless radius defined as the ratio of the radius 
at any point P to the disk rim radius, and where M is a centrifugal 
loading parameter defined by 


M = 


o + 
8 


V ? 4 
— nura- 


( 3 ) 


where v is Poisson’s ratio for the disk material. The primes in equa- 
tion (2) denote derivatives with respect to p. 
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v The boundary conditions used for all the solutions obtained herein 

are: 

* (a) The deflection at the shaft is zero. 

(b) The slope of the disk at the shaft is zero. 

(c) The total shear at the rim is zero. 

(d) 01he radial bending moment at the rim is zero. 

The equations expressing these boundary conditions as derived in appendix 
C are given by 


At 

P - 0 

II 

0 


At 

p = 8 

K 

II 

0 

\ 

At 

p = 1 

Y"i y" 

/ 

- (3 - v)Y' - (3 - v)Y = 0 

At 

p = 1 

Y" + vY* 

0 

11 

1 


These are by no means the only set of boundary conditions that could be 
used. For example, if a disk were to have radial rim loading, boundary 
* condition (d) which states that the radial bending moment at the rim is 

zero would be replaced by one that states that the radial stress at the 
rim is equal to the known rim load. 

The solution of this problem involves the solving of the fourth- 
order linear differential equation with variable coefficients (eq. (2)) 
subject to boundary conditions at two points (eq. (4)). For the case 
where M equals zero, the exact solution can be readily obtained. 

For M not equal to zero, recourse must be made to an approximate 
solution. The method used herein was to reduce the fourth-order differ- 
ential equation (2) to a set of four first-order differential equations 
and then to solve these by a Runge-Kutta finite -difference technique as 
described in appendix D. The method used for satisfying the boundary 
conditions is also described in appendix D. Once the differential equa- 
tion (2) is solved along with the chosen boundary-condition equations, 
the radial and tangential bending stresses o-^ and cr^ on the face of 
the disk due to the gyroscopic bending moments can be obtained by equa- 
tions (B20) as 

°R = K 


Orp — K 



* 
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where 

K = -12 mnft ( 2^ sin 6 + cos &) (6) 


KESUITS AND DISCUSSIOH 

Solutions were obtained using the finite-difference equations de- 
rived in appendix D for disks of constant thickness with the centrifu- 
gal loading parameter M varying from 0 to 50 and the ratio of shaft to 
disk radii varying from 0.1 to 0.3. In order to make the calculations 
for so many cases, a high-speed digital computer was used._ For the lim- 
iting case of M equal to zero, equation (2) reduces to the case of a 
stationary disk with normal loading. This problem is solved in closed 
form in reference 1 (p. 260). Figure 2 shows plots of dimensionless 
stress against dimensionless radius for both, the exact- and finite- 
difference solutions at this limiting condition of M equal to zero and 
for various values of (3 . It can be seen that excellent agreement has 
been achieved for all the cases investigated. 

Figure 3 shows plots of dimensionless stresses and deflections 
against dimensionless radii. Each of these plots is for one particular 
value of 6 and for a range of values of the centrifugal stress param- 
eter M. For a given set of operating conditions and disk geometry, the 
parameters C, M, and K are calculated from equations (l), (3), and 
(6), respectively. The stresses and deflections at any point in the disk 
can then be determined from the curves for the appropriate value of p . 

It should he noted from this series of figures that the dimensionless 
stresses and deflections both decrease as M increases . However, the 
actual stresses and deflections will always increase as M increases. 

For a given M, the effect of increasing p is to reduce the deflections 
and stresses as would be expected. These trends can also be seen in fig- 
ure 4, where the dimensionless stresses at the shaft and the dimension- 
less deflections at the rim of the disk are plotted against the centrifu- 
gal stress parameter M for various values of p. These plots are used 
to determine the most critical values of stresses and deflections for 
this type of problem. 

For an example of the vise of these figures in solving the problem 
of a parallel-sided disk, consider a disk with a thickness h of 0.50 
inch, an outside radius a of 9.0 inches, a shaft radius of 1.35 inches, 
an angular velocity co about its center of 2000 radians per second, an 
angular velocity ft and acceleration ft about a diameter of 1 radian 
per second and zero radians per second per second, respectively, and a 
mass dens ity for the disk material of 0.0008 lb-sec /in. , Poisson's 
ratio of 0.3, and a flexural rigidity of 0.345KL0 6 pound- inches . By 
vising these values to determine conditions along the radius 9 ** 0, the 



NACA UN 4=218 


5 


centrifugal stress parameter M (from. eq. (3)), the "bending stress param- 
eter K (from eq. (6)), and the deflection parameter C (from eq. (l)) 
are found to "be 12.6, 28 x 10 3 , and 0.275, respectively. For the value of 
M of 12.6 and the calculated value of p of 0.15, the dimensionless 
stresses at p = p and the dimensionless deflection at p = 1 are found 
from figure 4= to be 0.750, 0.225, and 0.580. When these dimensionless 
stresses are multiplied by the bending- stress parameter K, the radial 
and tangential stresses at the shaft are found to be 2000 and 6500 pounds 
per square inch, respectively. When the dimensionless deflection is mul- 
tiplied by the deflection parameter C, the maximum rim deflection is 
found to be 0.016 inch. 

It should be pointed out that these computed gyroscopic stresses are 
always vibratory stresses. As can be seen from the right side of equation 
(B17j, the gyroscopic load is a function of the angle 0 and, therefore, 
as the disk rotates, the load varies with the frequency o>/2it. The 
stresses, therefore, are vibratory stresses also with the frequency cs/2rt. 
The complete stress distribution for this disk is obtained by adding these 
vibratory bending stresses to the stresses due to the centrifugal loading. 

Although these calculations have been performed for the special case 
of a constant -thickness disk, the solutions for the more general equation 
(B19) for a variable-thickness disk may be performed in a similar way 
using the same finite-difference approach, the only difference being in 
the coefficients of the last of the four finite- difference equations (D8) 
and (D9) . The radial and tangential stresses R and T due to the cen- 
trifugal loading will have to be calculated in advance by any of the 
available methods such as in reference 2. It should also be noted that 
variations of Young's modulus E, due to either temperature or material 
variations through the disk, can be taken into account very easily by 
adjusting the flexural rigidity term D. 

In making calculations of this, type, the problem presents itself 
concerning what values of ft and ft are likely to occur during opera- 
tion. These values, in general, would be determined by the type of 
maneuvers the airplane or missile is capable of performing. However, 
the occurrence of air gusts and ground bumps may produce greater values 
of ft and ft. These effects are, of course, difficult to determine but 
some estimate would have to be made before a reasonably safe analysis 
could be made . 


CONCLUDING REMARKS 

Th e analysis of the gyroscopic loading on a constant -thickness disk 
has shown how this effect can be considered by the designer. For the 
case of a uniform-thickness disk, curves are presented that can be used 
for design purposes - The disk of variable thickness with temperature 
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gradient requires individual calculation, which can readily he performed 
by the method outlined in this report. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 19, 1957 
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APPENDIX A 
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SYMBOIS 

disk radius 

deflection parameter, 2 ntofla 5 -g sin 0 + cos 0 ^ 

flexural rigidity, Eh 3 /l2(l - v 2 ) 

Young’s modulus of elasticity 
disk thickness 

bending-stress parameter, - 12 nmfl s ^ n ® + cos 

centrifugal- loading parameter, ^ urn 2 a 4 ^ 

mass density 

normal load intensity 

radial stress due to centrifugal loading 

radial coordinate of point in disk 

tangential stress due to centrifugal loading 

disk deflection, function of p and 0 

dimensionless disk deflection, function of p coaly 

ratio of shaft to disk radius 

angular coordinate in plane of disk 

Poisson's ratio 

dimensionless radius, r/a 

radial bending stress on face of disk due to gyroscopic loading 

tangential bending stress on face of disk due to gyroscopic 
loading 


<p angular displacement of disk about a diameter x - x (fig. 5) 

Q or 9 angular velocity about disk diameter, rad/sec 

m or 0 angular velocity about disk center, rad/sec 
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APPENDIX B 


DERIVATION OP DIFFERENTIAL EQUATION FOR DEFLECTIONS OF VARIABLE- 
THICKNESS DISK SUBJECTED TO GYROSCOPIC LOADING 

If it is assumed that the load per unit area q acting on a plate 
is normal to its surface and that the deflections w are small in com- 
parison with the thickness h of the plate, the equation of equilibrium 
becomes (ref. 1, p. 85) 


^ - a Xy 

bx 2 by 2 &riy 


- q 


(Bl) 


where the bending moments M^., My, are given by 

Mjj = - ^(w^ + vWyy) 

My. = - B(Wyy 4" 

^ “ D(T - V) Vxy J 


(B2) 


where v is Poisson's ratio and D is the flexural rigidity of the 
plate. The bending stresses on the surface of the plate are 


6M* 




Cy - 9 

x h* 


a 

y h 2 


(B3 ) 


By substituting equations (B2) into equation (Bl) and letting D be a 
function of x and y, equation (Bl) becomes 

D(w + w + 2v ) + D (w„_ + vw )+D (w +vw ) + 
xxxx yyyy xxyy' xx s xx yy' yy v yy xx' 


2D xKxx + v xyy) + 2D y^yyy + v xxy } + 2(1 " v)l W = 5 (M) 


where the subscripts on D and w represent partial derivatives with 
respect to the subscript. If equation (B4) is converted to polar coor- 
dinates and it is assumed that D is only a function of the dimension- 
less radius p, equation (B4) becomes 
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\ D ' 
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vD" \ 

2 -n ) 

p D/ 


v ee 


D 1 vD^ 

P 2 D e D 

i 

+ ~4 w eeee 


v p + 



(B5) 


where the superscripts on D Indicate total derivatives with respect to 
p, and a is the outside radius of the disk. If equations (B2) are sub- 
stituted into (B3) and conversion to polar coordinates is made, the stress 


equations become 


°R “ 


= — 


6 D 


a 2 h 2 I PP 
6 D 

a 2 h 2 


[jpp + V (J w p + ^ w eej 

(| "p + 7 w ee + ™pp) 


\ 

> 


(B 6 ) 


The normal load intensity term q in equation (B5) may be analyzed 
as being composed of two separate loading terms . The first term, q^, is 

due to the so-called "gyroscopic forces," and the second term., q 2 , is 

due to the normal components of the stresses in the plane of the disk 
when the disk undergoes bending out of this plane (fig. l) . For the 
derivation of the q^ term, the coordinate system of figure 5 will be 
used. For an aircraft or missile making a turn, the turbine or compres- 
sor disk is rotating about its axis and at the same tine about a diam- 
eter fixed in space. Consider a point P on a disk as shown in figure 
5. The disk is rotating about its axis and also about a diameter that 
is at any time coincident with the X-axis . From reference to this fig- 
ure, it can be seen that the angular velocity cd is equal to 0, and the 
ang ula r velocity Q equals cp. Also, if it is noted from the figure 
that the magnitude of the vector r remains constant with time t, the 
orthonogal components of the acceleration of point P are given by 


Ajj. = (cos Q) 
dt 2 

fl 2 

A = — 5 (sin 9 sin cp) 

y dt 2 

rd 2 

A_ = — p (sin 6 cos cp) 


> 

J 


(B7) 
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where 0 is the angle "between the radius vector and the X-axis, and cp 
is the angle between the Z-axis and the projection of the radius vector 
cm the YZ-plane . The acceleration of point P normal to the plane of 
the disk is given by 


A n = aAjj. + bAy + cAg 


(B8) 


where a, b, and c are the direction cosines of the normal to the plane 
of the disk. For this case 


a = 0 
b = cos cp 
c = -sin cp 


1 


(B9) 


When equations (B7) and (B9) are substituted into equations (B8), the 
acceleration normal to the disk becomes 


Aj^ = r(’cp sin 9 + 2cp0 cos 0) 

The force per unit area on particle P is then given by 

qj_ = mhA^ = mhr(fl sin 9 + 2tx& cos 0) 


(BIO) 


(Bll) 


For the derivation of the second part of the normal load intensity 
q 2 due to components of the centrifugal forces when the disk undergoes 
bending, an element of the disk as shown in figure 6 will be used. Sum- 
ming up the Z-compcments of the forces on this element results in the 
following equation: 


h l 



SR 

3? 



Rrhc 


dS 3? 


+ T dr 


h / Sw 
r \§0 



mr dr d0 h — ^ = q 2 r dr d0 (B12) 

where R and T are the centrifugal radial__and tangential stresses, 
respectively, and the numerical subscripts on the thickness h and the 
partial derivatives of the deflection w represent their values at dif- 
ferent locations in the disk element as shown in figure 6. To relate 
both the thickness and the derivatives of the deflection at points 1, 

2, 3, and 4 to a common point at the center of the element, the first 
two terms of a Taylor series expansion is used, from which 


v 
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h + 


5h _dr 
2r 2 


. ah. dr 
*2 = h - 22“ 


> (B13) 


J 


where h 3 = h^ = h. Substituting equations (B13) back into equation (B12), 
dividing through by r dr d0 and dropping higher order terms, results in 


™ d 2 w ait Rh dh\ 2w Th a 2 w a 2 w /•- RnA \ 

42 i? a? ( 1 


The term. cS 2 w/dt 2 in equation (B14) can be written as 


bt 2 




2 a 2 w a 2 © Sw 

bd 2 Bt 2 


(B15) 


Substituting equation (B15) into equation (B14) and letting r = ap 
yields 


*2 


= ft W PP + ^2 (j** + T + h ' R ) V P + V 60 " ^“S© + W 0 ) (B16) 


When equations (B16) and (Bll) are added to give the total normal 
load, intensity term q, and this term is substituted into equation (B5) 
and terms involving like derivatives of w are collected, equation (B5) 
becomes 


w pppp + 2 


(l D l \ Tl , \ D' D" hRaf 

\p + D / PPP p 2 " ^ 2 + v ) p D " D D 


W PP + 


f-i .SL. , + M + h . B )| „ 

[p 3 p 2 D pD D \ p /] ( 


2(i -SL 

\p Z p 2 B 


v p©0 + 


_4 

■>P 


p + ^2 w ppse ■ 

3D' vD" hTa 2 2 
4 + P 2 D " p 2 D + D /" + 


• a%i 


— ^0000 + mi — = 2 wXl 


ha 5 ( h 


\2 


sin 0 + cos 0 




(B17) 
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If it is assumed that the disk is rotating about its center _ at a constant 
angular velocity during a change in direction of its axis, cd is equal to 
zero and a solution to equation (B17) may now be assumed to be of the 
form 


w « CY (B18) 

where Y is only a function of the dimensionless radius p, and where 

h 


C = 2ntoQ a 3 ^ B4n 0 + GOS 0 


•) 


Substituting equation (B18) into (B17) results in the total differential 
equation .. - .. 


Y"" + 2 




_3_ _ 3D_|_ vlT _ a^ 

IKR* + — + h'Rj 

p 3 p^) pD 15 

\ P /_ 


Y' - 


(: 


3 3D’ vD" hTa 2 2 ha 4 

— + — o — + nn> -p 


vp 4 p^ p 2 D p 2 D 
and the stress equations (B6) become 

3 - k r ' + ^' 


) Y - p 




T- -? Y ) 




where 


-12nn)fia 3 ( h „ . 

K = 5 Us 0 + cos 0 ) 


(B19) 


(B20) 


(B21) 
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For a so li d disk of constant thickness , the centrifugal radial and 
t ang ential stresses are given simply "by 

R = 3_+_v ^ 2 ( 3 . . p2) 

> (B22) 

J 

where the centrifugal stresses due to the disks rotation about its diam- 
eter have been neglected since 0 will always be much smaller than cd 
for all practical cases. All the derivatives of D and h with respect 
to p are equal to zero and under these conditions., equation (B19) re- 
duces to 


T 


3 + v 


ioi> 2 a 2 


( . 1 + 3v 9 | 

V- - 3T1T P v 



[\ + M(1 - p 2 )l Y" + 

& - M (rr^Jl 

\sr J 

p \ r 7 



(B23) 


where the centrifugal loading parameter M is given by 


Ms 


3 + v 2 4 h 
q wcn c & j 


(B24) 
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APPENDIX C 


f 


BOUNDARY- CONDITION EQUATIONS 


A set of four boundary conditions are needed for the solution of the 
fourth-order differential equation of the deflection of a disk under nor- 
mal loading. For a freely rotating single disk such as a turbine wheel, 
the boundary conditions will be taken as follows: 

(a) The deflection at the shaft is zero. 

(b) The slope of the disk at the shaft is zero. 

(c) The total shear at the rim is zero. 

(d) The radial bending moment at the rim is zero. 

The total shear V is given by (ref. 1, p. 259) 




(Cl) 


where 


a - d / . 1 d w 1 d 2 w\ 


and 


„ r-, ^ t> ( 1 S^w 1 dw\ 

Mpe - (l - v) d 3^ - ^2 Sej 

if equations (C2) are substituted into (Cl) and r = ap then 


v ~ / 1 d^w 1 3w 

a 5 W + p a P 2 ' P 2 


2 - v o w 




— d 2 wA 
Se 2 / 


(C2) 


(C3 ) 


If the total shear is equal to zero at the rim (or p = l) equation 
(C3) becomes 


b? w , S 2 w bw . \ c* 3 w 

— + - 3? + (2 - V) ^2 






( 3 - v) fg 

be* 


(C4) 


«r 
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The radial moment is equal to (ref. 1, p. 259) 


T> . / 1 dv 1 v\ 

** “ ' D [a? v V ^ P aP/J 


If r = ap and if this radial moment is equal to zero at the rim, 
tion (C5) becomes 

5*V 3w , ci^w 
0:= — + V ^ + V “ 






The last two boundary conditions are given simply by 


and 


w = 0 


Sw 


\ at p = p 


J 


(C5) 

equa- 

(C6) 


(C7) 


where 0 is the ratio of the shaft radius to the disk radius 
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APPENDIX D 




FINITE -DIFFERENCE EQUATIONS 

In order to obtain a solution to the general differential equation 
(B19), it is convenient to reduce this fourth- order equation to a set of 
four first-order equations as follows; let 

Y* - J 

Y" =. J' n K 

Y m =» K’ a L 

and equation (B19) becomes 

L' = E - AL - BK - CJ - DY (D2) 

where A, B, C, Q, and E axe the corresponding coefficients appearing 
in equation (B19). These equations can now be written in finite- 
difference form using a simple Eunge-Kutta method (ref. 3, p. 233). Con- 
sider a section of a disk as shown in figure 7 where there are n equally- 
spaced stations between p = (3 and p = 1. Consider the following Taylor 

series expansion for Y at the i + 1 station in terms of the value of 

Y at the i station 

A 2 

*1+1 - Y ± + ApY’ + + ••• (D3) 

If the first three terms of equation (D3) are used and it is assumed that 

the second derivative is approximated by 


(D4) 


If the first of equations (Dl) is substituted in (D5), equation (D5) 
becomes 


n 


Y* _ y* 
X i+1 X i 

Ap 


equation (D3) reduces to 


-i+l 


Y ± + 


Ap 

2 


C*i 


+ Y i+lJ 



(D5) 
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t 


^1+1 - + f C J i + Ji+i) 

and in a s imilar manner the values of J, K, and L at the 
i+1 station are obtained with the aid of equations (Dl) and 
(D2) 

J i + X - J i + f <-*i + "J+x) 

%+x ” *i + tt C 1 ! + L i+x) 

L i + x - L i-f B A i L i + B A + Vi + E i Y x> + 

C A i+l L i+l + B i-l-l K i+l ■+ °l+X J i+X + D i+X y i+x)] + 

X + E i+x) 


where the starred exponents represent first approximations to these val- 
ues. These terms will be obtained fl‘am the first two terms of the Taylor 
series or 



* 

%+! 


Y ± + ApJ ± 
J ± + ZspK ± 




LjL + i =* - Ap (Aj_L^ + J 


(D7) 


If the problem of a solid disk of constant thickness is to be solved, 
equation (B23) is substituted into the last of equations (D6) and (D7) and 
the values of Y, J, K, and L at the i+1 station can now be written in 
terms of the i station as 
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V 





T i+1 “ 

y i + ¥ ( J i + j i+i^ 




J l+1 “ 

J i + ¥ + ^+1^ 




Li+l “ 

% + ¥ fo. + Li+l^ 





-C.Hfc 

Li+l “ % 

**h 





p i+i L p i+i 

«» 

)j *1+1 

-tt-ter) 

¥ 

( Pi + 

p t+i^ 




y (m) 


“i+i 


^i+i 


-&HI 




**+i 


where the starred terms are given by 


L i+i - % + A> 


y* + i - y ± + 4% 

J i+i ■ J i + 



^1+1 “ % + A3l«i 



V (D9) 


C V s ) 

y i 

+ (p ± 5 

\ pi /. 




The solution of the differential equation (B23) has thus been re- 
duced to computation of a set of deflections and corresponding deriva- 
tives by means of the recurrance relations (D8) and (D9). It must be 
remembered that the solutions which satisfy equations (D8) and (D9) must 
also satisfy the boundary conditions (C8). Such a set of solutions can 
be obtained in the following manner: 

Let Y^, K^, and Y^, K^, *e two 

homogeneous solutions for Y, J, K, and L (that is, solutions for eqs. (D8) 
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co 

o 

C- 

■<* 


and (D9)) with the last terms in the L 1+1 and Lf +1 expressions, namely 
dp/2 (p^ + Pi+l) dp(p^), respectively, being equal to zero) and let 

(3) (3) (3) (3) 

y , J , K , and L be particular solutions for Y, J, K, and L. 
Furthermore, let these three solutions satisfy the boundary conditions 
that Y = Y' a 0 at p=p. This corresponds to Y^ and being 

zero. Since equation (B23) and the corresponding finite difference equa- 
tions (D8) are linear, the complete solution to equations (D8) can be 
written as 


aY? 1 - 1 + bY^ + y[ 3) 


_ ( 1 ) ._( 2 ) _( 3 ) 

(l) (2) (3) 

Kj_ » aK£ + bK^ + K> 


aL 


(i) 


+ bt 


? 2 > + 


l (3) 

X 




(DIO) 


J 


where a and b are determined so as to satisfy the two boundary condi- 
tions at pal. It should be noted that if the solutions were not chosen 
in order to satisfy two of the boundary conditions, two additional sets 

Y^ 5 ^, etc. would have to be included to satisfy the two 

boundary conditions at p = J3. To obtain a and b, equation (DIO) is 
substituted into the last two of equations (C8), or at pal. 


(al/ 1 ) + bl/ 2 ) + L^) + (aK^ + bK^ + K^) - (3 - v) RaJ^ + bj 1 - 2 ^ + J^ 3 ^) - (aY^ + bY^ + Y^jl 
\n n n/\n n n / a a n / \ n n n /J 


(aK^ 1 ^ + b^ 2 ) + B^ 3 ^) + v ["(ajt 1 ) + bj^ 2 ) + J^ 5 ^) - (a + bY^ + Y^J] 
\ n n n / [_\ n n n / \ n n n. /J 
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Tiie starting values for the three solutions that were assumed and that 
satisfy the boundary conditions at the first station p = 3 or i =» 1 
are 


id) = 

0 

yd) a 0 

id) _ 0 

jd) , 

0 

jd) a 0 

jd) a 0 

K w = 

0 

- 1 

*d) = 0 

t W. 

1 

a 0 

L d) a 0 


(D12) 


When these starting values are used in the recurrence relations (D8) and 
(D9) , values for Y, J , K, and L can be computed successively for all 
other stations. Once the values are computed for the last, or n station, 
the constants a and b are obtained from equations (Dll) and the com- 
plete solution is then given by equations (DIO). 


* * 
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Figure 1. - Loading on disk due to rotation about two 
axes. Arrow length, is proportional to gyroscopic 
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Dimensionless radius, p 
(a) Radial stress. 

Figure 2. - Comparison of exact dimensionless stress with stresses 
obtained by finite difference solution for various ratios of 
shaft to disk radii. Centrifugal load parameter, zero. 
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Dimensionless disk deflection, w/c Dimensionless tangential bending 
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Dimensionless radius, p 

(a) Ratio of shaft to disk radius, 0.100. 

Figure 3. - Variation of dimensionless radial stress, tangential 
stress, and deflection with dimensionless radius for various 
values of centrifugal loading parameter. 
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Dimensionless radius, p 

(b) Ratio of shaft to disk radius, 0.125. 

Figure 3. - Continued. Variation of dimensionless radial stress, 
tangential stress, and deflection with dimensionless radius for 
various values of centrifugal loading parameter. 
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(c) Ratio of shaft to disk radius, 0.150. 

Figure 3. - Continued. Variation of dimensionless radial stress, 
tangential stress, and deflection with dimensionless radius for 
various values of centrifugal loading parameter. 
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(d) Ratio of shaft to disk radius, 0.175. 

Figure 3. - Continued. Variation of dimensionless radial stress, 
tangential stress, and deflection with dimensionless radius for 
various values of centrifugal loading parameter. 


Dimensionless radial bending 
stress, a-ofv 



Dimensionless disk Dimensionless 

deflection, w/c tangential 


28 


NACA TN 4218 


J 




Dimensionless radius, p 

(e) Batio of shaft to disk radius, 0.200. 

Figure 3. - Continued. Variation of dimensionless radial stress, 
tangential stress, and deflection with dimensionless radius for 
various values of centrifugal loading parameter. 
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Dimensionless radius, p 

(f) Ratio of shaft to disk radius, 0.225. 

Figure 3. - Continued. Variation of dimensionless radial stress, 
tangential stress, and deflection with dimensionless radius for 
various values of centrifugal loading parameter. 
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Figure 3. - Continued. Variation of dimensionless radial stress, 
tangential stress, and deflection with dimensionless radius for 
various values of centrifugal loading parameter. 
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0 .20 .40 .60 .80 1.0 

Dimensionless radius, p 


(h) Ratio of shaft to disk radius, 0.300. 

Figure 3. - Concluded. Variation of dimensionless radial stress, 
tangential stress, and deflection with dimensionless radius for 
various values of centrifugal loading narameter. 
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0 10 20 30 40 50 


Centrifugal loading parameter, M 
(a) Radial stress. 

Figure 4. - Variation of dimensionless stress at shaft with 
centrifugal loading parameter for various values of ratios 
of shaft to disk radii. 








Dimensionless disk deflection, w/C 



0 10 20 30 40 

Centrifugal loading parameter 


(c) Tip deflection. 

Figure 4. - Concluded. Variation of dimensionless stress with 
centrifugal loading parameter for various values of ratios 
of shaft to disk radii. 
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